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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC).  The  Air  Force  project  manager  was 
Marshall  K.  Kingery.  The  results  were  obtained  by  ARO,  Inc.,  AEDC  Division  {a  Sverdrup 
Corporation  Company),  operating  contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force 
Station,  Tennessee.  Some  of  the  data  presented  in  this  report  were  obtained  under 
sponsorship  of  the  Army  BMDATC,  Huntsville,  Alabama,  for  Science  Applications,  Inc., 
El  Segundo,  California,  under  ARO  Project  No.  V41G-58.  The  work  was  done  under  ARO 
Project  Nos.  V32J-11,  V32S-P4,  and  V41G-59A.  The  manuscript  was  submitted  for 
publication  on  April  9,  1979. 

The  authors  wish  to  acknowledge  the  contribution  of  E.  O.  Marchand  formerly  of  ARO, 
Inc.  for  providing  calculations  of  real-gas  effects  presented  in  this  report.  These  calculations 
involved  adapting  the  method-of-characteristics  technique  of  Rakich  (Ref.  12)  to  blunt 
cones  for  both  real  and  perfect  gases. 
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1.0  INTRODUCTION 

The  possibility  of  real-gas  effects  on  the  aerodynamics  of  high-speed  vehicles  has  been  a 
concern  for  several  years.  To  investigate  such  effects,  a  series  of  shots  utilizing  blunt 
cones  was  made  in  the  von  Karman  Gas  Dynamics  Facility  (VKF),  AEDC,  1,000-ft 
Hyperballistic  Range  (G). 

Both  air  and  nitrogen  range  environments  were  used  in  the  investigation,  and,  as  it  was 
assumed  that  nitrogen  is  a  perfect  gas  for  the  present  test  conditions,  real-gas  effects  were 
deduced  from  differences  in  the  aerodynamic  measurements  made  in  these  two  test 
environments.  Measurements  were  desired  at  small  angles  of  attack  for  the  blunt  cones  that 
can  have  appreciable  force  and  moment  nonlinearities  even  at  the  small  angles.  Satisfactory 
data  reduction  of  the  small  amplitude,  nonlinear  flight  histories  required  a  reexamination  of 
data  reduction  procedures  used  in  free-flight  testing. 

The  purpose  of  this  report  is  to  present  the  aerodynamic  measurements  obtained  and  to 
describe  a  data  reduction  procedure  devised  to  reduce  the  nonlinear  free-flight  data. 

2.0  APPARATUS 


2.1  RANGE  G 

Range  G  consists  of  a  10-ft-diam,  1 ,000-ft-long  tank  that  is  contained  within  an 
underground  enclosure  (Fig.  1).  It  is  a  variable-density  aerodynamic  range  and  contains  53 
dual-plane  shadowgraph  stations.  Forty-three  stations  are  positioned  at  nominal  20-ft 
intervals,  yielding  an  840- ft  instrumented  length.  The  other  ten  stations  are  located 
approximately  10  ft  downrange  of  stations  5  through  10,  12,  13,  15,  and  16.  The  angular 
orientation  and  position  of  most  test  configurations  can  be  determined  to  within 
approximately  ±0.25  deg  and  ±0.002  ft,  respectively.  A  chronograph  system  measures 
intervals  of  flight  time  to  within  ±  2  x  10-7  sec.  The  range  vacuum  pumping  system  provided 
the  desired  nominal  range  pressures  of  20  and  50  torr  in  the  present  test  program.  The 
nominal  operating  temperature  of  the  range  is  76CF.  The  launcher  used  in  the  present  test 
program  is  a  two-stage,  light-gas  gun  with  a  2.5-in.-diam  launch  lube. 

2.2  CONES  AND  SABOTS 

The  test  configurations  were  10-deg  semiangle  cones  with  nominal  nose-to-base-radius- 
ratio  (l/d  values  of  0.07,  0.12,  and  0.3.  The  cones  tested  had  nominal  base  diameters  of  either 
1.0  or  1.5  in.  and  were  fabricated  of  aluminum  except  for  the  nosetips.  All  cones  had 
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nosetips  made  of  the  tantalum-tungsten  alloy  Fansteel  60T  to  minimize  possible  in-flight 
nose  geometry  changes  resulting  from  high  aerodynamic  heating  conditions. 

Sabots  used  were  of  a  conventional  four-component  design  and  were  fabricated  of 
Lexan®  .  All  cones  were  launched  in  an  uncanted  orientation  relative  to  the  sabot,  and  the 
initial  angular  disturbances  to  the  cones  were  those  arising  from  muzzle  effects  and  from  the 
cone-sabot  separation  process. 


A  e  d  c 
71-917 

Figure  1 .  Range  G. 


3.0  TEST  PROCEDURES 

The  initial  step  in  the  free-flight  range  testing  procedure  is  to  accelerate  the  model-sabot 
package  to  the  desired  velocity  with  use  of  a  launcher.  The  sabot  is  then  separated  from  the 
model  in  the  blast  tank  portion  of  the  range  by  aerodynamic  forces.  The  model  continues  on 
its  flight  through  the  range,  and  its  motion  histories  are  monitored  by  the  shadowgraph 
system.  Model  position-attitude-time  measurements  obtained  from  shadowgram  and 
chronograph  readings  are  used  as  inputs  to  the  different  data  reduction  computer  programs 
used  in  the  analysis  of  the  motion  histories.  In  this  test  program,  air  and  nitrogen  were  used 
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as  the  two  range  environment  fluids.  Shot  conditions  and  physical  characteristics  of  the 
cones  tested  are  tabulated  in  Table  1.  The  nominal  velocity  for  the  tests  was  16,000  fps. 
Nominal  Reynolds  number  levels  were  0.4  x  106  and  2.0  x  106. 

Laser  photography  stations  positioned  along  the  range  were  used  to  monitor  possible 
changes  in  the  nose  geometry  of  cones  arising  from  aerodynamic  heating  during  flight. 
Measurable  nose  geometry  changes  were  observed  only  in  the  downrange  portion  of  flights 
of  the  V'  =  0.12  cones  in  the  air  environment.  The  laser  photographs  permitted  detecting 
changes  in  the  effective  nose  diameter  of  a  model  to  within  about  .005  in.  Stability 
measurements  presented  for  these  cones  were  obtained  using  only  the  uprange  portion  of  the 
motion  histories  (to  shadowgraph  station  17,  the  first  320  ft  of  the  total  840-ft  instrumented 
length  of  the  range). 


Table  1 .  Test  Conditions  and  Physical  Measurements  for  the  Cones 


Shot 

Gas 

m 

gm 

d 

in. 

ly  *  104 

in. -lb-sec2 

eg 

V* 

fps 

Pr 

torr 

Tr 

deg  F 

1 

Air 

133.56 

1.491 

9.61 

0.6200 

0.1219 

16,645 

50  37 

71.1 

2 

Nitrogen 

130.96 

1.489 

9.38 

0.6201 

0.1209 

15,650 

49.14 

71.6 

3 

Air 

130.00 

1.488 

9.26 

0.6201 

0.1210 

15,334 

50.01 

71.6 

4 

Air 

130.83 

1.488 

9.48 

0.6200 

0.1210 

16,554 

50.25 

72.3 

5 

Nitrogen 

129.15 

1.487 

9.33 

0.6210 

0.1210 

15,899 

48.83 

72.5 

6 

Nitrogen 

130.70 

1.487 

9.40 

0.6201 

0.1210 

15,809 

48.53 

72.2 

7 

Nitrogen 

105.21 

1.489 

6.02 

0,5784 

0.3022 

15,764 

'48.70 

72.7 

8 

Air 

117.15 

1.488 

6.12 

0.5816 

0.3091 

16,001 

49.84 

72.2 

9 

Air 

105.41 

1.489 

6.07 

0.5741 

0.3088 

15,978 

50.15 

72.0 

10 

Nitrogen 

20.66 

0.9905 

0.531 

0.6303 

0.0707 

16,317 

19.56 

75.8 

11 

Nitrogen 

20.67 

0.9901 

0.534 

0.6304 

0.0706 

15,956 

19.54 

75.5 

Ref.  8 

Air 

21.50 

0.9928 

0.539 

0.6296 

0.0690 

15,187 

20.01 

74.0 

NOTE:  Roll  rates  of  cone  were  <  0.6  deg/ ft. 
“Midrange  Velocity 


4.0  DATA  REDUCTION 


The  numerical  integration  fitting  technique  (NIF  technique)  developed  in  recent  years 
(e.g.,  Refs.  1  and  2)  has  been  useful  in  reducing  free-flight  stability  data  for  axisymmetric 
vehicles  having  nonlinear  force  and  moment  characteristics.  Modifications  in  the  motion 
equations  discussed  in  Ref.  3  permitted  an  extension  of  the  NIF  technique  to  apply  to 
vehicles  with  asymmetries.  Such  asymmetries  include  the  primary  aerodynamic  asymmetries 
involving  appreciable  differences  between  CNa  and  CNiJ  and  between  |CniJ  and  |C,V|. 
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A  significant  feature  of  the  NJF  technique  is  that  motion  histories  of  more  than  one  flight 
can  be  fitted  simultaneously.  This  multiple  fit  feature  extends  the  capability  of  defining 
nonlinearities  by  permitting  the  use  of  different  flights  of  the  same  configuration  that  had 
experienced  large  differences  in  mean  amplitude  levels.  Experience  has  indicated  that  the 
amplitude  change  during  a  typical  range  model  flight  is  usually  not  sufficient  to  reduce  the 
nonlinear  aerodynamic  coefficients  adequately  from  a  single  flight.  This  is  particularly  true 
for  a  shot  having  a  small  initial  amplitude  of  about  5  deg  or  less. 

In  the  present  tests  of  blunt  cones  with  nonlinear  force  and  moment  characteristics, 
initial  cone  amplitudes  experienced  in  most  of  the  shots  were  4  deg  or  less;  hence,  a 
nonlinear  data  analysis  using  the  NIF  technique  was  therefore  not  feasible.  In  an  attempt  to 
provide  a  means  to  reduce  the  data  adequately,  quasi-linear  data  reduction  procedures  used 
in  the  past  were  reexamined,  focusing  on  reducing  center-of-pressure  data. 

For  an  ideal  linear  configuration,  the  moment  coefficient  can  be  written 

Cm  =  Cma  ‘  a  (1) 

where  CmQ  is  a  constant  for  different  values  of  a.  If  the  oscillatory  motion  histories  of  such  a 
linear  configuration  are  fitted  with  an  NIF  computer  program  (using  only  the  linear 
terms  of  the  program),  Cma  of  Eq.  (1)  is  one  of  the  parameters  defined.  The  amplitude  of 
the  motion  interval  fitted  is  of  no  concern  as  CmQ.  is  a  constant  and  Cm  can  be  defined  for  the 
a  values  of  interest. 

For  a  configuration  having  a  nonlinear  moment  that  can  be  described  with  a  cubic 
nonlinearity  (a  cubic  nonlinearity  is  normally  assumed  in  the  NIF  technique),  then 

Cm  —  Cm,a  •  ot  +  Cm3o  •  (2) 

If  a  motion  history  of  such  a  configuration  is  fitted  as  above,  using  only  the  linear  terms, 
then  the  motion-fitting  program  will  define  a  quasi-linear  moment  parameter  Cma.  The 
usefulness  of  Cmavalues  has  been  limited  in  the  past  because  of  lack  of  corresponding 
amplitude  information  available  from  the  motion  fit.  For  example,  if  measurements  of  Cma 
were  obtained  from  motion  fits  of  a  given  configuration  at  different  amplitude  levels  and 
they  were  approximately  constant,  then  such  measurements  would  have  confirmed  that  the 
moment  variation  for  the  configuration  was  approximately  linear.  However,  if  different 
values  were  obtained  for  Cma,  then  one  would  know  only  that  the  configuration  had  a 
nonlinear  moment  and  that  Cm,  which  is  of  primary  concern  in  a  nonlinear  system,  could 
not  be  defined  directly. 
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Murphy  (Ref.  4)  and  Rasmussen  and  Kirk  (Ref.  5)  presented  two  different  analytical 
methods,  each  devised  to  determine  nonlinear  coefficients  such  as  Cni]o  and  Cm3a  of  Eq.  (2) 
individually,  from  a  series  of  flight  histories  obtained  for  a  model  at  different  amplitude 
levels.  If  both  Clll|a  and  Cmia  could  be  evaluated  adequately,  then  the  nonlinear  variation  of 
Cn,  with  a  could  be  determined.  However,  as  is  discussed  later,  neither  method  would  be 
expected  to  be  applicable  for  the  small  amplitude  shots  of  the  present  investigation. 

If  a  CmQ  value  has  been  determined  from  a  fit  of  an  interval  of  motion  history  of  a 
vehicle  that  in  fact  has  a  nonlinear  moment  coefficient  defined  by  Eq.  (2),  then  it  follows 
that  there  is  some  effective  amplitude,  ae,  of  the  motion  interval  fitted  that  corresponds  to 
the  Cn,u  such  that  Cm  of  Eq.  (2)  can  be  w  ritten 

(Cm)a  =  ae  =  •  CKe  —  •  Q(e  +  Cmya  * 

The  obvious  problem  here  is  in  relating  ae  to  the  amplitude  of  the  motion  interval  fitted.  To 
accomplish  this,  the  ratio  (ac/am)  was  used  in  this  report.  Here,  am  is  the  average  of  the 
midrange  amplitude  of  the  fitted  interval  and  the  mean  of  the  initial  and  final  amplitudes  of 
the  fitted  interval.  If  ae  could  be  defined  adequately,  it  would  hold  particular  significance  in 
that  one  point  on  the  nonlinear  Cm  curve  could  be  defined  for  each  fitted  Cmf<  value  without 
the  problem  of  evaluating  both  Cm|a  and  Cm3ff  of  Eq.  (2). 

The  results  of  a  brief  study  of  the  ae  concept  for  realistic  free-flight  test  conditions 
indicate  that  the  ( ae/am )  ratio  tends  to  be  insensitive  to  the  magnitude  of  the  nonlinearity,  to 
the  type  of  motion  pattern,  and  to  the  amplitude  level  involved  for  a  fairly  wide  range  of 
test  conditions.  These  results  are  demonstrated  in  Figs.  2  through  4.  To  demonstrate  the 
characteristics  of  ae,  a  series  of  oscillatory  motion  histories  (representative  of  Range  G 
flights)  was  generated  for  different  levels  of  cubic  nonlinearities.  These  computer-generated 
“ideal”  motion  histories  (excluding  simulated  experimental  errors)  were  then  fitted,  using 
only  linear  terms  of  a  motion-fitting  program  to  obtain  corresponding  Cmo  values.  The  ae 
values  were  then  computed  using  the  relationship  obtained  from  Eq.  (3) 

ae  —  —  W 

Cm,  and  Cn,3a  are  both  knowm  parameters  (used  in  generating  the  motion  history),  and 
was  obtained  in  the  fitting  of  the  ideal  motion  history.  The  results  of  those  checks  are 
shown  in  Fig.  2.  The  plotted  moment  curves  represent  a  wide  range  of  cubic  nonlinearities 
and  correspond  to  the  generated  motion  histories  used.  A  calculated  moment  curve  for  one 
of  the  blunt  cones  of  the  present  investigation  is  also  shown  to  indicate  that  the  hypothetical 
moment  curves  used  here  are  reasonable.  All  determined  (ae/am)  values  fell  within  a  small 
band,  and  a  value  of  0.87  was  selected  as  a  representive  value  for  the  (ae/am)  ratio. 
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Figure  2.  Measured  («e/am|  ratio  for  different  levels 
of  cubic  nonlinearities. 

The  motion  histories  used  in  Fig.  2  (as  noted  in  the  figure)  correspond  to  a  fairly  narrow 
elliptic  pattern  with  an  ellipticity  (ratio  of  the  lengths  of  the  minor  axis  to  the  major  axis)  of 
0.25.  To  show  the  effect  of  motion  pattern  shape  on  the  (ae/am)  ratio,  a  series  of  additional 
motion  histories  was  generated  for  one  level  of  nonlinearity  for  a  band  of  ellipticity  values 
(representing  planar  to  circular  motion)  and  then  fitted,  again  using  only  linear  terms  of  the 
fitting  program.  Again,  ideally  generated  motion  histories  were  used.  The  results  of  those 
checks  are  shown  in  Fig.  3  in  which  (ae/otm)  is  platted  against  ellipticity.  The  measurements 
indicate  no  detectable  effect  of  ellipticity  on  the  (ore/am)  ratio  for  ellipticity  values  up  to 
about  0.5;  however,  (ae/am)  increases  with  further  increases  in  ellipticity  to  a  value  of  1.0 
for  circular  motion  as  would  be  expected.  Ellipticity  values  for  range  flights  are  generally  less 
than  0.3. 

A  demonstration  of  the  effectiveness  of  the  ac  concept  in  reducing  nonlinear  data  is 
shown  in  Fig.  4.  Motion  histories  were  generated  that  corresponded  to  three  different  types 
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Conditions  for  Motion  Histories 
Ge nerated  and  Fitted _ 


Figure  3.  Effect  of  motion  pattern  eHipticity  on  {ae/o;m). 

of  nonlinearities  (indicated  in  the  figure);  then  simulated  experimental  errors  (comparable  to 
errors  lhat  would  be  experienced  in  Range  G)  were  added  to  these  motion  histories.  These 
generated  histories  with  errors  were  then  Fitted  (over  different  amplitude  intervals), 
providing  a  Cn,Q  value  for  each  motion  fit.  The  corresponding  moment  curves  for  the  three 
types  of  nonlinearities  used  are  shown  in  Fig.  4.  Each  plotted  symbol  is  a  Cn1  value  obtained 
from  an  individual  motion  fit  using  the  relation  Cm  =  C  •  ae  and  0.87  for  (ac/am).  The 
small  deviations  noted  in  the  computed  Cm  values  indicate  that  Cm  can  be  determined  quite 
well  using  the  ae  concept  for  the  three  types  of  nonlinearities  and  different  amplitude  levels 
examined. 

The  results  of  the  above  investigation  indicate  that  Cm  values  can  be  obtained  quite  well 
using  the  ae  concept  from  low-amplitude,  nonlinear  free-flight  motion  histories,  given  the 
following  restrictions  that  (1)  the  test  configuration  is  statically  stable  at  zero  angle  of 
attack,  (2)  the  magnitude  of  the  nonlinear  Cm  curve  at  the  maximum  amplitude  of  the 
motion  interval  being  fitted  is  not  less  than  about  0.9  of  the  Cm  peak  value  of  the  interval,  (3) 
the  ellipticity  of  the  motion  pattern  is  less  than  about  0.5,  and  (4)  the  model  roll  rate  is  small. 
The  second  restriction  was  dictated  by  excessive  magnitude  of  deviations  obtained  in  the 
fitted  values  of  the  quasi-linear  force  parameter  (Cn^)  and  the  idamping  parameter  (Cni(.)  in 
motion  fits  of  cases  for  which  larger  decreases  in  the  moment  curve  existed.  The  magnitude 
of  nonlinearity  fora  moment  curve  with  an  increasing  slope  is  less  restrictive  in  motion  fits, 
and  the  ratio  (Cn,3a  •  cx^/ Cm|  )  can  be  at  least  three  times  larger  than  the  absolute  value  of 
this  ratio  for  the  moment  curve  with  a  decreasing  slope. 
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a,  deg 

Figure  4.  Capability  of  reducing  Cm  values  from  motion  histories 
using  (aB/am>  =0.87. 

Although  the  value  of  0.87  for  («e/am)  was  determined  using  cubic  nonlinearities, 
measurements  in  Fig.  4  indicate  that  the  insensitivity  of  this  parameter  permits  its  use  with 
noncubic  nonlinearities  similar  to  those  in  Fig.  4.  The  insensitivity  of  the  (ae/am)  ratio  under 
the  realistic  free-flight  test  conditions  of  the  present  study  is  not  too  surprising.  An  exact 
solution  was  obtained  by  Kirk  (Ref.  6)  for  a  system  having  a  corresponding  cubic 
nonlinearity  in  the  moment  for  the  basic  restraints  of  planar  motion  and  zero  damping. 
Although  the  method  of  Ref.  6  did  not  involve  the  use  of  the  ratio  (ae/am),  this  ratio  can  be 
calculated  from  presented  information.  Such  calculations  indicate  that  (ac/am)  is  essentially 
a  constant,  within  less  than  one  percent  (for  planar  motion  and  zero  damping)  for  the  above 
nonlinearity  restrictions. 

In  using  the  ae  concept,  the  normal-force  data  is  reduced  similarly  to  the  moment  data. 
Hence,  CN  =  CNo  «  ae,  where  CNcr  is  the  quasi-linear  normal-force  parameter  and  is 
determined  in  the  motion  fitting  process  similarly  to  Cma.  Previous  experience  in  reducing 
linear  free-flight  data  indicates  that  C1i1q,  a  function  of  the  frequency  of  the  oscillatory 
motion,  can  be  determined  very  well  (within  less  than  one  percent)  even  when  using  only  one 
cycle  or  so  of  the  motion  history.  Only  a  portion  of  a  typical  Range  G  flight  (which  usually 
has  three  to  four  cycles)  is  required  to  define  the  motion  frequency.  Hence,  the  larger 
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deviations  in  the  reduced  Cm  values  of  Fig.  4  indicate  the  limitation  in  defining  ae  for 
those  types  of  nonlinearities  used.  However,  this  is  not  the  case  with  free-flight  linear 
normal-force  data.  Here  CNq  is  obtained  from  the  fit  of  the  oscillatory  transverse  motion 
history  of  the  center  of  mass  of  the  vehicle,  and  basic  fitting  limitations  of  this  motion 
(transverse  motion  displacements  are  usually  small)  can  result  in  errors  up  to  ±4  percent  in 
CN(y  from  a  typical  range  flight.  Hence,  basic  motion  fitting  errors  can  be  expected  to 
predominate  over  any  ae-related  errors  in  determing  CN  values  from  use  of  the  quasi-linear 
normal-force  parameter  C^. 

Center-of-pressure  locations  can  be  determined  quite  well  from  free-flight  measure¬ 
ments  as  indicated  below. 

The  static  margin  (A)  of  a  vehicle  is  defined 

A  =  (Cm/CN)  =  (Cma/CNo)|inearsyslem  (5) 

and  the  center-of-pressure  location, 

cp  =  eg  +  A  (6) 

For  a  combined  error  in  the  (Cm/Crg)  ratio  of  5  percent,  a  3-percent  static  margin 
(referenced  to  body  length),  and  a  typical  eg  location  of  63  percent  of  vehicle  length,  the 
resulting  error  in  the  center-of-pressure  location  would  be  less  than  0.25  percent. 

This  method  [using  (ae/am)  =  0. 87]  should  be  adequate  for  handling  a  family  of  non- 
linearities  similar  to  those  shown  in  Fig.  4  and  is  not  restricted  to  just  a  cubic  nonlinearity.  It 
should  be  noted  that  the  method  is  applicable  in  reducing  free-flight  data  regardless  of 
whether  tests  are  made  in  aeroballistic  ranges  or  in  wind  tunnels. 

Although  the  ratio  (ae/am)  is  not  used  in  either  of  the  approximate  methods  of  Refs,  4 
and  5,  sufficient  information  is  presented  for  both  methods  to  permit  computing  some 
values  for  (ae/am).  Note  that  (ae/am)  values  corresponding  to  each  of  these  methods  for  a 
cubic  nonlinearity  and  planar  motion  are  in  agreement  with  the  results  of  the  current  study. 
However,  available  information  was  insufficient  for  such  a  comparison  for  nonplanar 
motion.  It  appears  that  the  linear  fitting  programs  that  have  been  involved  in  the  use  of  both 
of  these  methods  (Refs.  4  and  5)  incorporate  the  well-known  epicyclic  motion  equation 
defined  by  Nicolaides  in  Ref.  7.  In  that  type  of  a  motion  fitting  program,  the  quasi-linear 
parameter,  Cmc(,  obtained  in  a  fit  of  nonlinear  nonplanar  motion  can  be  dependent  on 
whether  the  model  roll  rate  is  measured  and  used  as  a  known  input  to  the  program  or  on 
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whether  the  roll  rate  is  permitted  to  float  without  restraint  in  the  fitting  process.  It  should  be 
noted  that  in  the  present  motion  fitting  investigation,  all  motion  fits  were  obtained  utilizing 
actual  roll  rates. 

The  drag  coefficient,  Cd,  was  obtained  from  the  relation 

CD  =  -(2m/eSV)  (dV/dx)  (7) 

where  dV/dx  is  the  mean  slope  of  the  curve  of  model  velocity  as  a  function  of  distance 
traveled,  and  V  is  the  mean  model  velocity  during  the  interval  of  the  flight  trajectory  being 
examined.  The  velocity  data  used  are  the  average  velocities  between  consecutive  stations  as 
obtained  from  the  measured  position-time  history  of  the  model.  The  term  dV/dx  is  obtained 
by  fitting  a  linear  equation  to  the  velocity-distance  data  by  means  of  a  least-squares 
procedure.  Equation  (1)  is  based  on  the  assumption  that  Cp  is  constant  during  the  flight 
interval  of  concern;  Eq.  (1)  has  been  found  to  provide  adequate  results  for  the  test 
conditions  of  the  current  program.  Since  the  angle  of  attack  of  a  model  will  normally  vary 
along  the  trajectory,  the  drag  coefficient  obtained  by  this  procedure  corresponds  to  the 
mean  total  angle  of  attack  experienced  by  the  model.  This  mean  angle  of  attack  is  defined  by 
the  relation 


d2  =  [1/(L2  -  L,)]  £«2dx  (8) 

where  L2  -  Lj  is  the  length  of  the  flight  interval  over  which  Cd  is  computed. 

The  zero  angle-of-attack  drag  coefficient  is  obtained  from  a  linear  extrapolation  of  a 
series  of  CD  measurements  (obtained  from  different  shots  made  at  similar  test  conditions)  as 
a  function  of  b2.  With  this  drag  reduction  procedure,  Cp  can  usually  be  reduced  by  using  the 
complete  flight  history  or  by  using  only  a  portion  of  the  flight.  This  procedure  is  useful  for 
examining  drag  measurements  in  which  possible  nose-blunting  effects,  resulting  from  high 
aerodynamic  heating  conditions,  may  be  involved. 

5.0  DATA  UNCERTAINTY 

The  estimated  total  uncertainties  (2 a  deviation)  associated  with  measurements  in  this  test 
program  are  as  follows: 


\ 
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Parameter  Uncertainty,  percent 


1.5 

Cm 

1.0 

CN 

4.0 

cp 

0.25 

+  Crtu 

10.0 

Any  systematic  errors  in  the  above  parameters  have  been  estimated  to  be  negligibly  small. 

6.0  RESULTS  AND  DISCUSSION 

Predictions  of  real-gas  effects  have  indicated  that  for  some  configurations  the 
aerodynamic  parameter  particularly  sensitive  to  such  effects  is  the  center  of  pressure.  For 
example,  predicted  center-of-pressure  locations  (cp)*  as  a  function  of  nose-to-base  bluntness 
ratio  (tA)  for  a  10-deg  semiangle  cone  are  shown  in  Fig.  5.  These  prediction*  were  obtained 
using  a  three-dimensional  method-of-characteristics  program  referred  to  herein  as  NASA- 
MOC.  For  A  values  of  about  0.12  (the  primary  A  value  of  the  present  investigation)  the 
predicted  curves  indicate  that  the  cp  in  an  air  environment  would  be  approximately  one 
percent  more  forward  than  the  cp  in  a  nitrogen  environment.  However,  a  reversal  occurs  at 
higher  A  values,  and  for  f  about  0.3  the  predicted  curves  indicate  that  cp  in  the  air  would  be 
about  one  percent  more  aft  than  the  cp  in  nitrogen.  A  shift  in  cp  of  this  magnitude  can  be  of 
much  concern  in  the  case  of  interceptor  and  maneuvering  vehicles  designed  with  small  static 
margins. 

All  of  the  aerodynamic  measurements  obtained  for  10-deg  semiangle  cones  are  presented 
in  Figs.  6  through  9.  Figure  6a  shows  both  the  cp  measurements  for  the  A  =  0.12  cone  and 
the  corresponding  calculated  cp  curves  (NASA-MOC)  for  nitrogen  and  air  environments. 
Each  experimental  point  shown  corresponds  to  a  given  range  shot  and  is  plotted  at  an  angle 
of  attack  equal  to  at  for  that  shot,  consistent  with  the  data  reduction  procedure  discussed  in 
Section  4.0.  The  cp  measurements  indicate  significant  real-gas  effects;  the  level  of  the  cp 
measurements  in  air  is  about  0.75  percent  more  forward  than  that  measured  in  nitrogen.  The 
calculated  curves  closely  predict  the  shift  in  cp  attributable  to  real-gas  effects;  however,  they 
predict  a  cp  about  0.5  percent  more  aft  than  the  corresponding  measured  cp  for  both 
nitrogen  and  air.  These  observed  differences  in  predicted  and  measured  levels  are  not 
understood  at  present.  It  should  be  noted  that  all  predictions  in  this  report  are  for  a  sharp 
edge  (zero  radius)  at  the  base  of  the  cones,  whereas  the  test  cones  were  fabricated  with  a 
nominal  edge  radius  of  0.03  in. 


“Provided  bv  E.  O.  Marehand,  formerly  of  ARO,  Inc. 
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Figure  5.  Predicted  center-of-pressure  location,  a  =  3  deg  (NASA-MOC). 

Because  of  the  significance  of  these  cp  measurements  relative  to  both  real-gas  effects 
and  differences  from  predicted  values,  the  range  measurements  were  scrutinized  for  possible 
sources  of  errors.  For  example,  the  sensitivity  of  the  center-of-pressure  location  with  nose 
bluntness  as  indicated  in  Fig.  5  would  at  first  appear  to  be  a  factor  of  concern;  however,  the 
nose  bluntness  was  measured  for  each  model  before  launch  (see  Table  1),  and  these 
measurements  indicated  that  the  model  noses  had  a  bluntness  ratio  slightly  greater  than  the 
design  specification  value  of  0.12.  Hence,  any  effects  of  bluntness  deviations  in  model 
fabrication  could  only  have  caused  decreases  in  the  differences  between  measured  and 
predicted  cp  values.  Another  possible  source  of  error  in  the  cp  measurements  is  the  nose¬ 
blunting  effecl  caused  by  aerodynamic  heating.  As  discussed  previously,  in-flight  nose 
blunting  was  observed  only  in  the  air  shots,  and  stability  data  were  reduced  using  only  the 
uprange  portion  of  each  of  those  three  flights.  Again,  any  undetected  nose  blunting  in  the 
uprange  flights  in  air  would  have  decreased  the  observed  real-gas  effects  measured  in  the  air 
and  nitrogen  shots.  The  results  of  such  examinations  indicate  that  any  errors  in  the  measured 
cp  levels  are  appreciably  smaller  than  the  differences  observed  in  Fig.  6a. 

The  cp  measurements  obtained  and  the  corresponding  predicted  values  are  shown  in  Fig. 
6b  for  the  V-  =  0.3  cone.  At  the  smaller  amplitudes  the  cp  measured  in  air  is  about  one 
percent  more  aft  than  that  measured  in  nitrogen.  The  calculated  curves  (NASA-MOC) 
predicted  the  real-gas  effects  well,  but  the  predicted  cp  levels  for  both  air  and  nitrogen  were 
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a.  \p  =  0.12 


b.  ^  =  0.3 

Figure  6.  Variation  of  cp  with  angle  of  attack. 

about  one  percent  more  forward  than  the  measured  values.  The  measured  cp  in  air  at  a  ~ 
11.5  deg  indicates  that  the  predicted  trend  of  aft  movement  of  the  cp  at  the  smaller 
amplitudes  reverses  above  5  deg.  This  is  supported  by  the  perfect-gas-calculated  cp  point  at 
ct  =  10  deg  in  Ref.  8.  Predictions  from  Ref.  8  (referred  to  as  NSVVC)  agree  reasonably  well 
with  the  NASA-MOC  perfect-gas  predictions  up  to  5  deg  and,  hence,  are  not  shown  herein. 
The  NASA-MOC  calculations  were  not  available  for  a  >  5  deg. 
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There  were  two  =  0.07  cone  shots  (nitrogen  environment)  made  primarily  to  assess 
possible  real-gas  effects  in  the  dynamic  stability  data  presented  in  Ref.  9;  these 
measurements  are  discussed  later  in  the  report.  Although  these  cones  experienced  large 
amplitude  motion,  cp  measurements  were  obtained  for  them  and  are  shown  in  Fig.  6c  along 
with  a  cp  data  point  for  a  i/-  =  0.07  cone  (Ref.  10).  The  calculated  cp  curve  (air  environment) 
shown  in  Fig.  6c  agrees  with  the  data  from  the  low-amplitude  shot  discussed  in  Ref.  10;  a 
small  real-gas  effect  is  indicated  for  the  \j/  =  0.07  cone  by  the  air-nitrogen  predictions.  The 
cp  values  measured  for  the  two  high-amplitude  shots  indicate  that  the  cp  moves  forward 
slightly  with  increasing  angle  of  attack.  The  perfect-gas  cp  prediction  (NSWC)  at  a  =  10  deg 
is  in  reasonable  agreement  with  the  high-amplitude  experimental  points. 


c  =  0.07 
Figure  6.  Concluded. 

The  Cm  and  CN  data  for  the  =  0.07,  0.12,  and  0.3  cones  are  shown  in  Fig.  7. 
Significant  real-gas  effects  are  apparent  in  the  moment  data  for  both  the  ^  =0.12  and  0.3 
cones  at  small  angles.  No  appreciable  real-gas  effects  were  observed  in  the  normal-force  data 
obtained.  The  ^  =  0.12  moment  data  show  appreciable  nonlinear  Cm  variations  measured 
for  both  air  and  nitrogen  environments.  The  high-amplitude  data  points  (a  >  5  deg)  are  in 
reasonable  agreement  with  the  small-amplitude  data.  The  calculated  curves  predict  the  shifts 
in  the  force  and  moment  data  arising  from  real-gas  effects  well;  however,  some  of  the 
predicted  levels  are  noticeably  displaced  from  the  measured  Cm  levels.  No  normal-force 
data  were  obtained  for  the  f  =  0.12  cone  shots  because  of  the  cones’  experiencing  small 
transverse  motions,  resulting  primarily  because  the  cones  had  larger  static  margins  than  are 
desirable  for  good  normal-force  measurements.  The  cp  measurements  for  these  shots  were 
reduced  using  predicted  Cm  values  at  the  corresponding  «c  values.  This  is  a  sensible  approach 
since  normal-force  predictions  are  generally  quite  good  for  cones.  (The  agreement  of 
predicted  and  measured  CN  values  for  the  \p  =  0.07  and  0.3  cones  at  small  amplitudes  should 
be  noted.) 
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3.  \f/  -  0.07 

Figure  7.  Moment  and  normal-force  coefficients  as  functions 
of  angle  of  attack. 

Results  of  an  aeroballistic  range  investigation  made  several  years  ago  of  the  aerodynamic 
characteristics  of  10-deg  semiangle  cones  at  high  velocities  were  published  in  Ref.  9. 
Measurements  of  free-flight  damping-in-pitch  derivatives  therein  presented  indicated  an 
appreciable  increase  in  damping  with  increasing  Mach  number  between  M  =  8  and  16.  The 
free-flight  measurements  were  made  in  an  air  environment;  there  has  been  some  concern 
that  real-gas  effects  were  reflected  in  those  measurements.  Figure  23  of  Ref.  9  comparing 
the  damping  measurements  has  been  replotted  in  Fig.  8.  Also  shown  in  Fig.  8  are  two 
damping  measurements  obtained  in  the  present  investigation  for  the  two  ^  =  0.07  cones. 
Those  shots,  made  in  a  nitrogen  environment,  agree  with  the  previously  obtained  free-flight 
measurements  made  in  air  and,  hence,  indicate  no  significant  real-gas  effects  in  the  damping 
measurements.  After  the  two  ^  =  0.07  cone  shots  were  made  in  nitrogen,  Ericsson  and 
Reding  (Ref.  1 1)  showed  that  the  increasing  damping  trend  with  Mach  number  between 
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M  =  8  and  16  is  realistic  when  viscous  crossflow  effects  are  accounted  for  in  the  predicted 
damping  values.  It  should  be  noted  that  the  tests  of  Ref.  9  covered  a  wide  range  of 
amplitudes  and  indicated  no  detectable  nonlinearities  in  the  damping  measurements  for  a  < 
0.07. 


Tunnel-Free  Oscillation  (Sting-Supported  Models) 
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Figure  8.  Comparison  of  damping  derivatives  (Fig.  23  of  Ref.  91. 


Drag  measurements  for  the  three  blunt-cone  configurations  are  shown  in  Fig.  9.  No 
detectable  real-gas  effects  on  the  drag  coefficient  were  observed.  Of  additional  interest  are 
the  shots  for  the  \p  =  0.12  cone  made  in  an  air  environment.  As  noted  previously,  these  shots 
experienced  nose  blunting  (caused  by  aerodynamic  heating)  during  their  nights.  The  effect 
of  nose  geometry  changes  resulting  from  aerodynamic  heating  and  corresponding  ablation 
effects  on  CD  apparently  tends  to  be  compensating  (see  Fig.  9b,  drag  measurements  reduced 
over  total  length  of  flights)  for  the  present  test  conditions.  Such  compensating  effects  on  CD 
have  been  observed  in  other  test  programs. 
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In  summary,  (he  test  measurements  indicate  significat  real-gas  effects  in  the  current 
center-of-pressure  and  moment  data.  Consistency  in  these  measurements  indicates  that  the 
devised  procedure  for  reducing  small-amplitude,  nonlinear,  static  stability  data  is  quite 
feasible.  Because  of  the  large  inertia  values  involved  in  the  present  tests  of  the  ^  =  0.12  and 
0.3  cones,  obtaining  damping  derivatives  for  these  cones  was  infeasible. 


52,  deg2 

a.  i  =  0.07 


62,  deg2 


b.  !/■  =  0.12 


5 2  deg2 

c.  \j/  =  0.3 


Figure  9.  Drag  coefficients  for  the  blunted  cones. 
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7.0  CONCLUDING  REMARKS 

A  free-flight  investigation  of  real-gas  effects  on  the  aerodynamics  of  blunted  10-deg 
semiangle  cones  was  made.  The  nominal  test  velocity  was  16,000  fps,  and  the  nominal  Ref 
levels  were  0.4  x  106  and  2.0  x  106.  The  test  models  were  -  0.07,  0.12,  and  0.3  cones. 

Real-gas  effects  as  evaluated  from  measurements  using  different  test  environments  (air 
and  nitrogen)  were  appreciable  in  the  center-of-pressure  and  moment  data.  For  example, 
shifts  in  the  cp  position  of  up  to  one  percent  of  the  cone  length  were  measured.  No 
significant  real-gas  effects  were  observed  in  the  normal-force  or  drag  data.  A  comparison 
using  damping  data  previously  published  indicated  no  apparent  real-gas  effects  on  the 
damping  deiivatives  for  the  ^  =  0.07  cone.  Measured  changes  in  the  center-of-pressure  and 
moment  data  arising  from  real-gas  effects  were  predicted  well  with  the  NASA-MOC 
method. 

The  aerodynamic  measurements  made  in  the  present  test  program  indicate  that  the  new 
procedure  devised  to  reduce  nonlinear  static  stability  data  from  small-amplitude  motion 
histories  is  feasible.  The  data  reduction  procedure  is  applicable  for  free-flight  tests  made  in 
either  aeroballistic  ranges  or  wind  tunnels. 
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NOMENCLATURE 


Cd 

Cm 

CN 


CftMqd/iV)  +  C|r-(rtd/2V) 

Cn/J 

c~ 

ma 

Cnna,  Cjn  ja 


CNa»  CN|j 


eg 


cp 


d 

ly 

t 

m 

Po 


Pr 


Rep 


Mean  drag  coefficient  for  flight 
Moment  coefficient 
Normal  force  coefficient 
Damping  in  pitch  derivatives 
Moment  derivatives 
Quasi-linear  moment  derivative 
Nonlinear  moment  coefficients 
Normal  force  derivatives 
Quasi-linear  normal  force  derivative 

Center  of  gravity,  measured  from  nose  of  the  cone  and  referenced 
to  the  cone  length 

Center-of-pressure  location  measured  from  the  nose  of  the  cone 
and  referenced  to  the  cone  length 

Cone  base  diameter 

Moment  of  inertia  (relative  to  a  transverse  axis) 

Cone  length 
Cone  mass 
Initial  roll  rate 
Range  pressure 

Reynolds  number  based  on  free-stream  conditions  and  cone  length 
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Tr 

V 


a 


«e 


am 


6 


^o> 

'/'o 


CjJ 

SUBSCRIPT 

x.xx? 

<30 


Range  temperature 
Cone  velocity" 

Total  angle  of  attack  (small  angles)  or  sine  of  total  angle-of-attack 
Effective  angle  of  attack,  see  Eq.  (3) 

Mean  amplitude  of  motion  interval  being  fitted 

Total  angle  of  attack 

Mean  total  angle  of  attack,  see  Eq.  (8) 

Euler  angles  (initial  values) 

Time  derivatives  of  the  Euler  angles  (initial  values) 

Cone  nose  bluntness  ratio  (nose  radius/base  radius) 

Model  frequency 


Moment  reference  point  location 
Ambient  conditions 

NOTE:  All  coefficients  and  derivatives  are  referenced  to  the  cone 
base  area.  The  reference  length  for  all  moment  coefficients  and 
moment  derivatives  is  the  cone  length,  except  where  indicated 
otherwise. 
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